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Abstract 

Vacuum expectation values for one Wegner- Wilson loop representing a moving 
quark-antiquark pair are calculated in four-dimensional Euclidean and Minkowski space- 
time. The calculation uses gluon field strength correlators with perturbative gluon 
exchange and non-perturbative correlations from the stochastic vacuum model. The 
expectation value of a Wegner- Wilson loop forming a hyperbolic angle in Minkowski 
space-time is connected by an analytical continuation to the expectation value of the 
Wegner- Wilson loop in Euclidean four-space. The obtained result shows how confine- 
ment enters into the light-cone Hamiltonian for valence quarks independently of the 
chosen model. 
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1 Introduction 



One of the challenges in quantum chromodynamics (QCD) is the relativistic bound state 
problem. In the light-cone Hamiltonian approach [1] light-cone wave functions can be con- 
structed in a boost invariant way. It is necessary to have reliable light-cone wave func- 
tions if one wants to calculate high energy scattering, especially exclusive reactions. Many 
parametrizations assume separability of the dependence on the longitudinal momentum frac- 
tion and transverse momentum which is very unlikely since the two momenta are coupled 
in the kinetic energy operator. Various approaches have been tried to compute such wave 
functions. One can use the usual equal time Hamiltonian [2] and transform the resulting 
wave functions into light-cone form with the help of kincmatical on-shell equations. The 
light-cone Hamiltonian in a string picture is formulated in ref. [3]. More ambitious is the 
construction of an effective Hamiltonian including the gauge degrees of freedom explicitly 
and then solving the bound state problem. For mesons this approach [4, 5] still needs many 
parameters to be fixed. Attempts have been made to solve the valence quark wave function 
for mesons in a simple Hamiltonian with a two-body potential [6]. 

A necessary input is an adequate potential for the light-cone Hamiltonian. For the equal 
time Hamiltonian and heavy quarks the calculation of Wcgner- Wilson loops gives the form of 
the non-perturbative potential for long distances. The correlator model [7] allows to calcu- 
late vacuum expectation values of gauge invariant Wegner- Wilson loops using perturbative 
and non-perturbative field strength correlation functions as input. One computes the loop 
expectation value (Wr[C]) in terms of a gauge invariant bilocal gluon field strength correlator 
integrated over minimal surfaces by using non-Abclian Stokes' theorem. Then the matrix 
cumulant expansion in the Gaussian approximation is applied. 

The basic object of the correlator model is the gauge invariant bilocal gluon field strength cor- 
relator Ffj^^pf^{Xi,X2,0;Cxj^o,Cx2o)- The strings C^io, C^^o connect the coordinates Xi, X2 
in the correlation function of the two-field strengths to a common reference point O. We 
define 
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The gluon field correlator has a perturbative (P) and a non-perturbative {NP) component. 
The stochastic vacuum model is used for the non-perturbative low frequency background field 

and the perturbative gluon exchange for the additional high frequency contributions. The 
most general form of the correlator respecting translational, Lorentz and parity invariance 
reads in Euclidean space [7] 
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with G2 = (^■^G'^p,y{0)G'^^^,{0)^ as the gluon condensate. The term proportional to n is the 
non-Abehan confining part F^^"^ of the correlator, in contrast, the tensor structure Fj^J^^'^ is 
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characteristic for Abelian gauge theories and does not lead to confinement. The correlation 
functions are a simple exponential of range a: 



The calculation of a Wegner- Wilson loop along the imaginary time directions gives the heavy 
quark- ant iquark potential with color-Coulomb behavior for small and confining linear rise 
for large source separations [7]. Since the computation of the VEV for one Wegner- Wilson 
loop can be done completely analytically, also other orientations of the loop can be chosen, 
e.g. a loop where the quark-antiquark pair moves along the z-direction. By transforming 
to Minkowski space-time the dependence of the interaction potential on longitudinal and 
transverse separation of the pair can be obtained this way. In section 2 we describe the 
calculation in Euclidean space-time, in section 3 in Minkowski space time and in section 4 
we derive the potential in a light-cone Hamiltonian for valence quarks. 

2 Vacuum expectation value for a tilted Wegner- Wilson 
loop in Euclidean space-time 

The vacuum expectation value (VEV) of a tilted Wegner- Wilson loop represents a moving 
quark-antiquark pair 



The group generators are in the fundamental representation of SU{3), g is the strong 
coupling constant and V symbolizes path ordering of the closed path C in space-time. The 
loop C with spatial extension Rq and temporal extension T has the following parametrization 
in four- dimensional Euclidean space-time (Fig. ^ 
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{ut^, ue[-T/2,T/2]} 
{T/2t^ + vRor^, ve[0,l]} 
{-ut,, + Ror^, ue[-T/2,T/2]} 
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and the parametrization of the surface 



= ut^ + vRor^, u e [-T/2,T/2]- v G [0, 1], 
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Figure 1: Configuration of the Wegner- Wilson loop in Euclidean space-time. 

The expression for the Wegner- Wilson loop simplifies with the help of the Casimir operator 
in the fundamental representation C2(3) = = 4/3 



(iy[C])G = exp 



C2(3) 



-Xs 



where Xss is the double area integral of the correlation function over the surface 
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The lengthy calculation of Xss is standard and follows the fines of reference [7] , we will give 
here only these parts which are relevant to understand the calculation of the tilted loops. 
One gets for the non-perturbative confinement component 
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j dvi j dv2 j dui j du2D{Z^). (11) 
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Correlated points on the surface have the distance Z = {ui — U2)t + {vi — V2)Ror. The 
geometry of the loop orientation enters via the factor a 



o? = \ — cos^ 0sin^ Q. 
The confining Xss has the following final form: 
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At large distances Roa » 2a one recognizes that the confining interaction leads to a VEV 
of the tilted Wilson loop which is consistent with the area law Ro 



<W[C]> = e-'"^""^ (14) 
= ' (15) 
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where a is the string tension [7] and the area is obtained from 



Area ^ TRo du dv^ i ] i -[-7^-7^] 16 
J-1/2 Jo \ \ du J \ dv J \ du dv J 

= TRoa. (17) 

The non-confining Xss functions give the short range attractive quark- antiquark interaction 
from massive correlator and gluon exchange 

- "t™ 3(iV2-l) ^H^J ^^^^ 

^ _ V^exp(-mci?oa) 

T^oo AirRoa ^ ' 

In the limit of straight loops {9 = 0) all results agree with previous calculations [7]. 

3 Vacuum expectation value for one Wegner- Wilson 
loop in Minkowskian space-time near the light-cone 

In this section the vacuum expectation value (VEV) of one Wcgncr- Wilson loop near the 
light-cone is computed in Minkowskian space-time. As before wc use the correlator model 
for the non-perturbative low frequency background field and perturbative gluon exchange for 
the additional high frequency contribution. The path of the color dipole in four- dimensional 
Euclidean space-time is represented by a light-like QCD Wegner- Wilson loop: To accom- 
plish the transition from Euclidean to Minkowski space-time we have to make the following 
replacements: 



Xi — > ix° 



(20) 



Here X — (X, X4) is the Euclidean space-time point and x — (x^, x) the Minkowskian vector. 
Of course, also the Euclidean correlation functions of the field strengths have to be changed 
to the Minkowskian correlation functions [8]. 

The loop C with spatial extent Rq and temporal extension T is now placed in four- dimensional 
Minkowski space-time. The quark-antiquark pair is moving with velocity /? 

/3=^ (21) 



5 





/ / 


/ 


D / 







Figure 2: Configuration of tlie Wegner- Wilson loop in Minkowski space-time, 
and the hyperbolic angle ip defines the boost (Fig. |2I). 

The loop has the same parametrization as before but with Minkowskian vectors r'^ and 
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For the Wegner- Wilson loop in Minkowski space-time we define Xss in the same way as in 
ref. [8] 
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The phase factor Xss is given as the double area integral over the surface 
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Xss ■■= -i^ Jda^^ixi) J daf"'{x2) F^^p^{xi,X2,0;C^,o,C^^o)- (24) 
s s 

The expression for the Wegner- Wilson loop is the same as before eq. (jH)) and eq. Q. We 
use the above convention in order to be consistent with the notation of ref. [8]. In the course 
of the calculation we find 



X 



NPc 



12(iV, 



- I da^^^ix,) I daP''{x2) iD^z^la") {g^^g,, - g^^g,,} , (25) 



where the changed metric tensor in Minkowski space comes from replacing the Euclidean 
metric tensor by the Minkowski tensor in the correlation function. Using the Minkowskian 
vectors r, t with the following properties 
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= cosh^ ^ — sinh^ ip = 1 

= r^r^ = — 1 

{t-rf = {tf.r'') ■ (fr^) = cos^ sinh^ ^ 

■ r"^ - it ■ rf = -(1 + cos^0sinh^^) 



one obtains 



T-^oo i 

The three-dimensional correlation function D^^\p^) in Minkowski space is obtained by ana- 
lytical continuation of the Euclidean function. Since the argument of the function is given 
by the magnitude of a three-vector the analytic continuation leads to the same modified 
Bessel function which was obtained in Euclidean space, cf . also Appendix B of ref . [8] . The 
resulting Xss is real. The geometry enters via the factor 

alj^l + cos^ sinh^ ijj (26) 

which is consistent with the analytical continuation of the Euclidean expression a — 1 — 
cos^ 0sin^ 9 into Minkowski space by transforming the angle 9 ^ iip. This analytical contin- 
uation is similar to the analytical continuation used in high energy scattering [11-13] where 
the angle between two Wilson loops transforms in the same way. 
The different contributions to Xss read: 



Roa 



M 
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Xss - -^1™ 3(Ar2_i) ■ RoO^M ■ K2 [—^ ) (27) 
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We remark that the same calculation with a time-like vector t-t — and a quark separation 
r'^ oriented orthogonally to t, i.e. r ■ t = gives the well-known result that the expectation 
value of the loop equals unity as has been shown in the original references on high energy 
scattering. In the formulation of high energy scattering [8] with the correlation function 
summed to all orders in the S-matrix the respective phase factors from single loops totally 
cancel. So the here obtained expectation value does not change the results of the previous 
calculations even if one goes away from the light-cone. 
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4 Using the VEV of the Hght-Hke Wilson loop to de- 
rive the quark-antiquark potential in the light-cone 
Hamiltonian 

The exponent giving the expectation value of the Wilson loop acquires a new meaning for 
a tilted loop. In order to interprete the result of the preceding section one must define the 
four- velocity of the particles described by the tilted loop 

u^^ {^,0^,^(3). (29) 
With the help of the four- velocity we can rewrite the loop as: 

The hue integral of the gauge potential acts as a phase factor on a Dirac wave function t/j 
which splits up into a leading dynamical component and a dependent component i/j^. 
For very fast quarks the mass term and transverse momenta are negligible compared with 
the energy and longitudinal momentum. In this eikonal approximation the Dirac equation 
of the leading component decouples from the small component: 



zd_i,+ = P-,(A-)^+ (31) 
= 9A-^+. (32) 

With /9 ~ 1 the phase factor in the tilted Wilson loop integrates A~ and leads to a VEV for 
the loop containing P~„t = " P^)\pot 

< Wr[C] >= e-^^(^°-^')l-*^. (33) 

The light-cone potential energy arising from the confining part of the correlation function 
has a term of order 0(p+) = O(^), where P+ = -^{P^ + P^) is the light-cone momentum 

P^ot = ^ (aRo^cos{<l>r + sin(0)772) . (34) 

Terms involving transverse momenta and masses of the same order 0(-p+) are not included 
in the loop as it has been calculated. Two of these terms give the standard kinetic energy 
term of free particles, which contributes to the total light-cone energy. Terms with spin 
cannot be obtained from the straight Wilson loop and are not discussed here. We introduce 
the relative + momentum and transverse momentum k±^ for the quarks with mass /i. By 
adding the above "potential" term to the kinetic term of relative motion of the two particles 
we complete our approximate derivation of the light-cone energy P" 

To derive the light-cone Hamiltonian we multiply P~ with the plus component of the light- 
cone momentum P+ = 73^"^° ~^ ^^^^ P'^/M = \/27 to eliminate the boost 
variable from the Hamiltonian. Further we follow the notation of reference [14] and introduce 
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the fraction ^ = k'^ / with |^| < 1/2 and its conjugate the scaled longitudinal space 
coordinate \/2p = P^xs as dynamical variables. For our configuration the relative time of 
the quark and antiquark is zero 

+ kl) 



Ml = 2P+P- = jj^:^ + 2a^p2 + M^xl (36) 

We have obtained the light-cone Hamiltonian from the confining interaction in a Lorentz 
invariant manner, because the variables ^, p, k±_ and x±_ are invariant under boosts. The 
valence quark light-cone Hamiltonian has a simple confining potential. The magnitude of 
the confining potential is set by the string tension a. The effective "distance" of the quarks 
is given by scale-free light-cone longitudinal distance and the transverse distance multiplied 
by the bound state mass. The above equation agrees in the limit of one-dimensional motion 
with the equation for the yo-yo string derived in ref . [14] . If there is only transverse motion 
(,^ = 0), the confinement has the usual form which is seen by using M ^ 2fi. To solve the 
operator one can go over to M^. Minimizing with respect to s one can replace the 
square root operator. Final self consistency must be reached with a guessed mass eigenvalue 
Mo 



^ {^^' + kl) ^ 1 L y + Mlxl 



+ s\. (37) 



1/4-^2 ' 2 V s 

The other non-confining potentials from the abelian correlator and the perturbative gluon 
exchange can be worked out similarly and one gets for the complete valence Hamiltonian 

,.2 + w (2g^M^e-^\ ^(7(l-«;)r2^ r , 

1/4-^2 ' I 47rr / kMott ^Ma' ^ ' 

with the dimensionless variable 



r = yJp'^ + M'^xl. (39) 

The relative weight of non-perturbativc non-abelian and abelian contribution is fixed by 
K = 0.7 in the paramctrization of the correlation function [8]. We have used a to parametrize 
also the abelian non-confining potential, instead of giving the full expression with the gluon 
condensate. Of course, the abehan part of the potential does not confine. The best way 
to find the two-body wave function is to keep ^ in the momentum representation and the 
transverse direction in configuration space x±. It is to be expected that in this approximation 
the pion mass is not described correctly. Firstly, the spin structure of the meson is not 
refiected in the spin independent expression above. Secondly, one expects quark self energy 
corrections [15] . The typical mass scale of a meson estimated from a trial solution with a wave 
function ^(,^, x^) = A cos (^7r)e~^i/^o comes out to be 1 GeV. The problem of the pion has to 
be addressed separately, since on the light-cone the mechanism of chiral symmetry breaking 
needs special care and is of particular interest. In the approach above confinement plays an 
important role in contrast to Nambu-Jona-Lasinio effective models which give an adequate 
description of spontaneous chiral symmetry breaking but do not include confinement. 
The confining interaction in the light-cone Hamiltonian was derived in the specific model of 
the stochastic vacuum. But it also can be inferred from the simple Lorentz transformation 
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properties of the phase in the Wilson loop and a lattice determination of the tilted Wilson 
expectation values. In this respect the final Hamiltonian is model independent. 
The inclusion of confining forces in the initial and final state wave functions can put all 
scattering cross sections calculated with the stochastic vacuum model on a much safer base 
since wave functions and cross section are derived consistently. For low the long distance 
part of the photon wave function matters strongly and confinement is important cf. [16]. 
Especially the diffractive cross section has a sizeable contribution from large dipole sizes and 
a correct behaviour can only be expected when the problem of the large dipole wave function 
is treated adequately. A very useful extension of the above calculation is the coupling of 
the initial qq state to higher Fock states qqg with gluons or fragmentation of the original 
qq state into qqqq states which can also be estimated from Wilson loops near the light-cone 
in Minkowski space. On the problem of fragmentation the Lund model [17] has been very 
successful and it is interesting to see how the above model calculation fares in comparison. 



Acknowledgments: We thank A. Hebecker for bringing this problem again to our attention 
and our colleague H.G. Dosch for a helpful discussion. 

References 

[1] S. J. Brodsky, H. C. Pauh and S. S. Pinsky, Phys. Kept. 301 (1998) 299 
iarXiv:hep-ph/9705477j . 



[2] 
[3] 
[4 

[5 
[6 

[7 
[8 
[9 



S. Simula, Phys. Rev. C 66 (2002) 035201 |arXiv:nucl-th/0204015| . 



V. L. Morgunov, V. I. Shevchenko and Y. A. Simonov, Phys. Lett. B 416 (1998) 433. 
M. Burkardt and S. Dalley, Prog. Part. Nucl. Phys. 48 (2002) 317 



| arXrv:hep-ph/01 120071 . 
S. Dalley and B. van de Sande, Phys. Rev. D 67 (2003) 114507 |arXiv:hep-ph/0212086 



T. Frederico, H. C. Pauh and S. G. Zhou, Phys. Rev. D 66 (2002) 116011 
|arXiv:hep-ph/0210234 . 



A.I. Shoshi, F.D. Steffen, H.G. Dosch and H.J. Pirner, Phys. Rev. D 68 (2003) 074004. 
A.I. Shoshi, F.D. Steffen and H.J. Pirner, Nucl. Phys. A 709 (2002) 131-183. 



O. Nachtmann, High Energy Collisions and non- perturbative QCD, \hep-ph/960936^ 
September 1996. 



[10] Andreas Kramer, Nichtperturbative Effekte in elastischer Hadron-Hadron-Streuung, 
Heidelberg University, Diss., 1991. 

[11] E. Meggiolaro, Z. Phys. C 76 (1997) 523 )ar Xiv:hep-t h/9602104l. 

[12] E. Meggiolaro, Nucl. Phys. B 625 (2002) 312 |arXiv:hep-ph7ori0069j . 



10 



[13] A. Hebecker, E. Meggiolaro and O. Nachtmann, Nucl. Phys. B 571 (2000) 26 
|arXiv:hep-ph/9909381|. 

[14] W. A. Bardeen, 1. Bars, A. J. Hanson and R. D. Peccei, Phys. Rev. D 13 (1976) 2364. 

[15] Y. A. Simonov, Phys. Lett. B 515 (2001) 137 |ar Xiv:hep-ph/0105141| . 

[1 6] H. G. D osc h, T. Gousset and H. J. Pirner, Phys. Rev. D 57 (1998) 1666 
|arXiv:hep-p'h/9707264| . 

[17] B. Andersson, Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol. 7 (1997) 1. 



11 



